Stretched exponential behavior in remanent lattice striction 
of a (La,Pr)! 2 Sr! 8 Mn 2 07 bilayer manganite single crystal 
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We have investigated the time dependence of remanent magnetostriction in a 
(La,Pr)i.2Sri.gMn207 single crystal, in order to examine the slow dynamics of lattice distor- 
tion in bilayered manganites. A competition between double exchange and Jahn- Teller type 
orbital-lattice interactions results in the observed lattice profile following a stretched exponential 
function. This finding suggests that spatial growth of the local lattice distortions coupled with 
e 9 -electron orbital strongly correlates with the appearance of the field-induced CMR effect. 

PACS numbers: 75.47.Lx,75.50.Lk 



The discovery of colossal magnetoresistance (CMR) ef- 
fect in doped manganites with perovskite structure has 
stimulated considerable interest for the understanding 
of their physical properties 1 . Though the insulator to 
metal (IM) transition and its associated CMR are well ex- 
plained on the basis of the double exchange (DE) model, 
it is pointed out that the dynamic John- Teller (JT) ef- 
fect due to the strong electron-phonon interaction, plays 
a significant role in the appearance of CMR as well as 
the DE interaction 2 ^. Furthermore, Dagotto et al pro- 
pose a phase separation model where the ferromagnetic 
(FM) metallic and antiferromagnctic (AFM) insulating 
clusters coexist and their model strongly supports recent 
experimental studies on the physics of manganites 4 . 

The bilayer manganite Lai.2Sri.sMn207 exhibits a 
paramagnetic insulator (PMI) to ferromagnetic metal 
(FMM) transition around T c =^120K and its associ- 
ated CMR effect. In comparison with cubic mangan- 
ites, the MR effect of the compound under considera- 
tion, due to its layered structure, is enhanced by two 
orders of magnitude, at 8T, around T c . It is well 
known that Pr-substitution on the La-site leading to 
(Lai_£,Pr z )i.2Sri. 8 Mn 2 C>7 causes an elongation of the c 
axis length in contrast with a shrinkage of the a(b) axis, 
resulting in a change of the e ff -electron occupation from 
the d x 2_ y 2 to d 3z 2_ r 2 orbital 6 - 7 . For z—0.6, a strong 
suppression of the e g -electron bandwidth leads to a PM 
insulating ground stated. However, the FM metastablc 
phase with a metallic conduction is easily obtained by the 
application of a magnetic field. The field induced first- 
order IM transition is further accompanied by a mem- 
ory effect of magnetoresistance, magnetization, magne- 
tostriction and magnetothermal conduct ivity£*2ii£. This 
memory effect may be explained by a schematic picture of 
the free energy exhibiting two local minima correspond- 
ing to the FMM and PMI states, as shown in the inset 
of Fig. 1(a). In the absence of magnetic field, the sys- 
tem stays in the stable PI state but, upon field appli- 
cation, the FMM state becomes metastable against the 



PMI state. After removing the field, the system should 
go back to a pure PMI ground state through the various 
kinds of mixed states after a very long time. Recently, 
Gordon et aim have reported that the time dependence of 
remanent magnetization of (La,Pr)i.2Sri.sMn207 follows 
a stretched exponential function, which is closely related 
to a magnetic frustration between FM double exchange 
and AFM super exchange interactions. In several reports 
concerning cubic manganites, the magnetization and re- 
sistivity relaxation effects have been discussed within the 
framework of the phase segregation between FM metal 
and AFM charge ordered insulating statesii*i2*i2ii4. In 
addition, there are a few reports on magnetic and trans- 
port properties in doped bilayer manganites which extend 
our understanding of a spin-glass like phase related to the 
phase separation model-' 1 - . However, to our knowledge, 
slow dynamics of a lattice relaxation in doped mangan- 
ites has not been reported so far. 

In this paper, we report the time dependence of re- 
manent magnetostriction in the Pr-substituted bilayer 
manganites. A frustration of lattice deformation between 
double exchange and Jahn- Teller type orbital-lattice in- 
teractions plays a crucial role in the observed slow relax- 
ation going from the FM metal to PM insulator. 

Single crystals of (Lao.4,Pro.6)i.2Sri. 8 Mn 2 07 were 
grown by the floating zone method using a mirror fur- 
nace. The calculated lattice parameters were reported 
earlier in 8 . The dimensions of the z=0.6 sample are 
3.4x3 mm 2 in the a6-plane and 1mm along the c-axis. 
Magnetostriction, both in the a6-plane and along the o 
axis , was measured by means of a conventional strain 
gauge method at the Tsukuba Magnet Laboratory, the 
National Institute for Materials Science and at the High 
Field Laboratory for Superconducting Materials, Insti- 
tute for Materials Research, Tohoku University. The 
magnetization measurements were made by using a su- 
perconducting quantum interference device magnetome- 
ter at Iwate University. A magnetic phase diagram in 
the (H, T) plane was established from the magnetic mea- 
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surements on the z=0.6 crystal as shown in the inset of 
Figl.(b)Siil. The phase diagram is separated into three 
regions labeled as the PMI, FMM and mixed phases (the 
hatched area in the (H,T) plane). In the mixed phase, 
the PMI and FMM clusters coexist, which is character- 
ized by hysteresis in the magnetization curves. Thus, we 
tried to look for a relaxation effect for the z=0.6 crystal 
at temperatures selected between 20 and 30K, covering a 
large to small hysterises region. After cooling the sample 
down to the respective temperatures in a zero field, the 
field was applied along the c-axis up to 5T at a sweep rate 
of 0.2 T/min and was then decreased down to zero at the 
same rate. A dotted horizontal line in the phase diagram 
represents such an experimental procedure crossing the 
phase boundary separating the hysteresis and FM areas. 
Finally, magnetostriction data were recorded as a func- 
tion of time just after the field was switched off. 
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FIG. 1: Magnetostriction both in the afo-plane and along the 
c-axis, for single crystalline (Lao.4,Pro.6)i.2Sri.gMn207. A 
schematic picture of the free energy with two local minima 
corresponding to the FMM and PMI states is shown in the 
inset of Fig. 1(a). The inset of Fig. 1(b) represents a magnetic 
phase diagram in the (H,T) plane established from the mag- 
netic measurements carried out on the z=0.6 crystal. 

Magnetostriction data both in the ab-plane and along 
the c-axis , for the (Lao.4,Pro.6)i.2Sri. 8 Mn 2 07 single crys- 
tal are shown in Fig. 1 as a function of magnetic field ap- 
plied parallel to the c-axis at selected temperatures. At 
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FIG. 2: The remanent c-axis normalized magnetostriction 
of single crystalline (Lao.4,Pro.6)i.2Sri.gMn207 as a function 
of time, just after the field was switched off. The solid 
lines correspond to a fit by a stretched exponential function 
, dL c (t) / dL c (to)=l-exp[— (t/r) 13 ], where r and j3 represent 
the characteristic relaxation time and exponent, respectively. 
(a)r=6.9 x 10 4 s and /3=0.35 at 20K,(b)-r=6.0 x 10 3 s and 
/3=0.42 at 25K and (c)r=1.7 x 10 2 s and /3=0.48 at 30K. 



20K, the value of dL c (H) / L C (Q) exhibits an abrupt de- 
crease around 3T , while the value of dL a b(H) / L a b(0) 
rises at the same field value. Here, the value of dLi(H) 
is defined as Li(H) — Lj(0). In other words, the o 
axis length suddenly shrinks around 3T in contrast with 
the expansion of the a(b) axis accompanied by the field- 
induced IM transition. Combining dL a \,{H) / 'L a b(0) and 
dL c (H)/L c (0), a volume shrinkage associated with the 
IM transition is deduced, indicating a drastic suppres- 
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sion of the local lattice distortion associated with the 
transition from the localized to itinerant stated. The re- 
sultant volume striction remains invariant down to zero 
field although its magnitude depends on the respective 
temperatures. In the insulating state of hole doped bi- 
layer manganites, there exists a lattice disorder between 
the Mn 3+ C>6 and Mn 4+ 06 octahedral, with and with- 
out the JT local lattice distortion, respectively^. As we 
pass T c from the insulating state, the Mn-0 bond length 
disorder is increasingly suppressed and, in the metallic 
state, it becomes less inhomogeneous than that found in 
the insulating state. Giant magnetostriction of the z=0.6 
crystal is understood on the basis of a strong suppression 
of the local lattice disorder due to the field-induced FM 
phase, where the double exchange interaction is favor- 
able against the orbital-lattice interaction. Moreover, the 
polarized neutron diffraction study on the z=0.6 crystal 
suggests that the field-induced metallic state, where the 
field is parallel to the c-axis, is characterized by a high 
population of the d3 Z 2_ r 2 orbital of Mn 3+ — . This finding 
is consistent with the c-axis distortion based on the local 
JT effect in the insulating state, preferring occupation of 
the d^ z 2_ r 2 orbital through the orbital-lattice coupling. 

Next, let us show in Fig. 2 the remanent c-axis mag- 
netostriction data of the (Lao.4,Pro.6)i.2Sri.gMn207 sin- 
gle crystal as a function of time, just after the field 
was switched off . Here, the value of dL c (t)/dL c (to) 
denotes the normalized dL c (t) = [L c (t) — L c (0)} by 
di c (to) = [L c (t < 0) - L c (0)] where L c (t < 0) and 
L c (t = 0)correspond to the virgin and initial values, be- 
fore the application of field and just after the removal of 
field, respectively. If the system fully relaxes to the PM 
virgin state , dL c (t) should approach an equilibrium value 
of dL c (to). It should be noted that, at higher tempera- 
tures, above 30K, the remanent striction instantaneously 
decays just after removing field, while at low tempera- 
tures (5K), the value of dL c (t) remains constant for 1 
day. At 30K, the lattice striction rapidly rises after re- 
moval of the field, then it relaxes by a half of the total 
value up to the relaxation time (=1.7 x 10 2 s) and finally 
restores the ground-state value. Furthermore, at lower 
temperatures, it takes a much longer time for dL c to de- 
cay to the ground-state value, as shown in Fig. 2. The 
temporal profile of the remanent magnetostriction fol- 
lows a stretched exponential function dL c (t) / dL c (to)=l- 
exp[— (t/r) 13 ], where r and (3 represent the characteristic 
relaxation time and exponent, respectively. The fitted 
parameters, r and (3, are plotted as a function of tem- 
perature in Fig. 3(a). Here, the temperature dependence 
of r is well described by a thermally activated function 
t = ,4exp(A/fcT) where A = 2.0x 1CT 2 and A = 30meV. 
Moreover, the value of the exponent is roughly fitted by a 
linear function of the form (3 = 3 x 10~ 2 + T/66.7. Such 
a deviation from the simple exponential decay needs a 
further treatment on the basis of a multiple relaxation 
time approximation to describe the observed relaxation 
effect. The transition probability (r _1 ) from the FMM 
metastable to the PMI stable states is proportional to 



exp(-5U / kT) , where SU represents potential barriers be- 
tween the metastable state and the local maximum. A 
wide distribution of the potential barriers is related to a 
stretched exponential decay, leading to our understand- 
ing of the phase separation between FM and PM clusters. 
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FIG. 3: (a) remanent magnetostriction. Characteristic relax- 
ation time and exponent, r and (3, respectively, plotted as a 
function of temperature. The former is well fitted by a ther- 
mally activated function, while the latter follows a T-linear 
function. For comparison, in the inset of Fig. 3(a), the fitted 
parameters to remanent magnetization data are also shown, 
(b) c-axis normalized magnetization of the same crystal as a 
function of time after removing the field. r=3.1 x 10 6 s and 
/3=0.16 at 20K, r=1.2 x 10 4 s and (3=0.22 at 30K. In the in- 
set of Fig. 3(b), the data collected for M C (H) before measuring 
the remanent magnetization, are also shown. 



For comparison, the c-axis magnetization of the same 
crystal is presented in Fig. 3(b) as a function of time af- 
ter removing the field. In the inset of Fig. 3(b), are also 
shown the data collected for M C (H) before measuring the 
remanent magnetization. In the inset of Fig. 3 (a), the pa- 
rameters, t and (3, obtained by fitting the magnetization 
data, are plotted as a function of temperature 22 . In a 
similar way, the activation energy of the relaxation time 
is estimated to be 33 meV and the exponent is fitted to 
the T-linear function/? = 4 x ICT 2 + T/167 . We notice 
that the characteristic r of L c is about by two orders of 
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magnitude smaller than that of M c on the same crystal, 
accompanied by the higher exponent . 

Here, we comment some differences in the relaxation 
processes between the magnetostriction and magnetiza- 
tion although both decay curves follow a similar stretched 
exponential function. First, in the view point of a frus- 
tration among the competing interactions, both the sit- 
uations are different. The remanent magnetostriction 
slowly relaxes because of the competition between the 
double exchange and JT type lattice-orbital interactions. 
The former interaction causes a suppression in the lo- 
cal lattice distortion along the c-axis through the itiner- 
ant state, but the latter favors the lattice deformation 
through the local John- Teller effect. On the other hand, 
in a previous study carried out on the z—Q.Q crystal, the 
temperature variation of the magnetization exhibited a 
steep decrease in the ZFC scan below T S =28K, which is 
probably associated with a spin-glass stated. It is thus 
expected from the spin-glass like behavior of M(T) that 
the decay process of remanent magnetization is closely 
related to a magnetic frustration between FM double 
exchange and AFM super exchange interactions. Fur- 
thermore, we point out the discrepancy in the fitted re- 
laxation time between L(t) and M(t) curves as shown 
in Fig. 3. A gradual drop in M(H) curve at 20K with 
decreasing field is quite contrast with field-independent 
behaviors of L(H) and p(H) at low temperatures during 
the field sweep to zeroS. The similarity in remanent val- 
ues between L(H) and p(H) data strongly suggests that 
the system still remains a metallic state with a volume 
shrinkage just after removing the field. As one of reasons 
for the difference between M(H) and L(H) (or p(H)), 
we think that the nucleation of magnetic domains in the 
metallic region probably gives a decrease in M(H) at low 
fields. 

Recently, Y.Tokunaga et al. have revealed 
the existence of magnetic domain structures in 
(Lao.4,Pro.6)i.2Sri.8Mn207 single crystal using a 
magneto-optical imaging technique 2 ^. It is true that 



the formation of domain walls never produce a much 
longer time of magnetic relaxation than the lattice 
relaxation. However, some of domain walls and/or each 
spins are pinned at lattice defect sites, giving a long 
lifetime of magnetic domains. Another scenario for the 
disagreement is concerning the existence of magnetic 
polarons in the insulating state, causing a magnetic 
contribution. In particular, in the hole doped bilayer 
manganites, Zhou et al., proposed that Zener polarons 
occupying two Mn sites forms a ferromagnetic cluster 
in the paramagnetic regime above T c , to explain the 
pressure effect on thermoelectric power—. In addition, 
it should be noted that the Debye- Waller (DW) factor 
along the Mn-apical oxygen bond in the insulating state 
is greater— than the DW in the metallic one, at most 
by a factor of 2. In other words, the difference in the 
local lattice inhomogeneity between the metallic and 
insulating states is not as large as that in M between 
the ferromagnetic and paramagnetic states. Finally, we 
make some comments on the relaxation phenomena of 
remanent resistivity in manganites. The resistive relax- 
ation studies support that the insulator-metal transition 
and its associated CMR effect are well described on 
the view point of both the phase separation and its 
related percolation behavior of metallic phase 11,25 . The 
spatial growth of the local lattice distortions coupled 
with e ff -electron orbital strongly correlates with the 
appearance of the field-induced CMR effect through the 
spatial evolution of percolating paths. 

In summary, we have reported the time dependence of 
remanent magnetostriction in the Pr-substituted bilayer 
manganites. A frustration between double exchange and 
Jahn- Teller type orbital-lattice interactions plays a cru- 
cial role on the observed slow relaxation of lattice going 
from the field-induced FMM to PMI states. 
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